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ABSTRACT 

1— It is shown that the upper bound for a in the general solutions of spherically sym- 



\^J ' metric vacuum field equations (gr-qc/9812081, A = 0) is nearly 10^. This has been 



Q"^ ■ obtained by comparing the theoretical prediction for bending of light and precession 

^r . of perihelia with observation. For a significant range of possible values of a{a > 2) 

y^ , the metric is free of coordinate singularity. 2— It is checked that the singularity in the 

^D I non-static spherically symmetric solution of Einstein field equations with A (JHEP 

04 (1999) Oil, a = 0) at the origin is intrinsic. 3— Using the techniques of these 
two works, a general class of non-static solutions is presented. They are smooth and 



O 

O^' finite everywhere and have an extension larger than the Schwarzschild metric. 4— 

$H . The geodesic equations of a freely falling material particle for the general case are 

rr| solved which reveal a Schwarzschild-de Sitter type potential field. 

>■ 

X 



PACS numbers: 04.20.Jb, 04.70.Bw, 98.80 Hw. 



I. INTRODUCTION 

According to Birkhoff' s theorem the metric for a vacuum spherically symmetric gravita- 
tional field with A = 0, has a unique Schwarzschild form 

ds' = (1 _ ^)rft2 _ (1 _ ^)-irfr2 - r\dd^ + sin^ Odct?) (1) 

and with A 7^ it has a unique Schwarzschild-de Sitter form 

rfs2 ^ (1 _ ^ _ ^r^)de -(I- — - -rY'dr^ - r^dO' + sin^ Odcj)^) (2) 

r 3 r 3 

Eq. ([l|) has a coordinate type singularity at r = 2M and an intrinsic singularity at r = 
while Eq. (^ has two coordinate singularities at r ^ 2M and r ^ (f)^? and an intrinsic 
singularity at r = 0. The intrinsic singularity is irremovable and this is indicated by diverging 
the Riemann tensor scalar invariant M 



48M 



2 



7^'^ bcdTZa '''' = —^ (3) 



For case A = 0, a general class of solutions has been obtained which has the form p[ 

ds^ = (1 - ^^ )de - (1 - ^^ )-Mr^ - (r + aM)2(rf^2 ^ ^^^2 ^^^2) ^4) 

r + aM r + aM 

where a is an arbitrary constant. In Sec. II we primarily discuss the apparent objection that 
they are subspaces of the Schwarzschild metric. Though it seems only a linear change of 
variables, but a does affect the curvature of space-time. Then it is shown that staticness 
is not a necessary initial condition if ggg is independent of time. Subsequently by deriving 
the equation of precession of perihelia and bending of light in a gravitational field and 
comparing them with observational measurements, we conclude that a may take small as 
well as large values up to 10^. Even though this may solve the singularity problem, but recent 
observations of type la supernovae indicate the existence of a positive cosmological constant 
0]. From other side, it has been shown that in the presence of cosmological constant, using 
static reference coordinate, is not suitable 0]. 

For case A 7^ 0, a non-static solution for this system has been proposed which has the 
following form ||^ 



ds' = {^ — '- — '-- — ^ — ^ — '- — '--)de 



-eVt* 



1 - ^ - f p2)2 + Mp2 + (1 _ 2M _ y^) 



+ r\de^ + sin^ 9dcj)'^) 



(5) 



where p = re v 3*. The result which evidently is free of any singularity for r 7^ 0, is singular 
at r = 0. In Sec. Ill we show that it is indeed an intrinsic singularity. In Sec. IV, by making 



use of the presented techniques in p,|[ , a general class of non-static solutions will be obtained 
which has the form 



1 _ 2M _ A 2)2 + 4A 2 ^^ _ 2M _ A 2) 



-eVt* 



;i_Mi_V)2 + ^^2+(l_2M_A^2) 



2 

+ {r + aMf{de^ + sin^ edcf)^) 



'dr^ 

(6) 



where here p = (r + aM)e^ ^^ . These solutions are smooth and finite everywhere even at 
r = 0. We show that they have an extension larger than the Schwarzschild metric. Obvi- 
ously they should be checked for completeness before we may call them non-singular. Sec.V 
deals with solving the geodesic equations for a freely falling material particle in the general 
case and results a potential field which though is very large aX r ^ but it is finite. Finally, 
derivations of Eq.(|2^), Eq.(|26|) and Eq. (p9D which being used in Sec. II, are presented in 
Appendixes A and B. 



II. CASE A = , a / 

Since Eq. (^ transforms to Eq. (|lD by simply replacing r' = r + aM with the range 
of r' > aM, this may cause a confusion that (|) is a subspace of ([^). The proof of 
completeness usually for a pseudo-Riemannian manifold is not an easy task. We will show 
the flaw in this argument by a Riemann counter-example. Taking R^ as a two-space of all 
points with coordinate r, 9 such that the metric is 

ds" = dr^ + r^de^ (7) 

where 6* = is identified with 9 = 2t[ and the point r = is included. This plane is complete 
and non-singular. Also R'"^ is a two-space of all points with coordinate r, such that the 
metric is 

ds^ = dr^ + (r + afde^ (8) 

where the range of r and 6 is the same as R^. If we transform r' = r + a Eq. (|]) gives 

ds^ = dr'^ + r'^de^ , and a < r' (9) 

which apparently this means R'^ C R^. We show that indeed this is not valid. Let's con- 
sider a subspace of R^ and R'^ by restricting r < p. The surface area of R'^{r < p) is Tip^ 
while the surface area of R''^{r < p) is vr(p^ -|- 2ap). This means for finite p we always have 
R'^{r < p) C R'^{r < p). If we take the limit p — > cxd then we get _R^(r < p) -^ R^ and 
R''^{r < p) ^ R'^. Since R^ is complete and from other side R'^{r < p) G R''^{r < p), 
there is no way except to conclude that R^ = R'^. This counter-example shows that how 
the conclusion that the general solutions are subspaces of the Schwarzschild metric may be 



impulsory. Indeed the Schwarzschild solution Eq.([l|), is a special case of the general solutions 
Eq.(^ for the case a = 0. We must notice that both the Schwarzschild and the general solu- 
tions are in the same Schwarzschild coordinate, which manifestly have different forms. Any 
transformation to the new radial coordinate r' = r + aM requires the Schwarzschild metric 
be written in this new coordinate too, which means we should replace r by r' — aM. Thus 
in this new coordinate they will have different forms too. As we know the Schwarzschild 
metric is an exterior solution which is only valid for r > 2M and incomplete. This is true 
for the general solutions but r > (2 — a)M gives an extension larger than the Schwarzschild 
solution and even with a > 2 the range of validity is r > 0. In the following we try to show 
how the spaces of general solutions for a > 2 are complete. 

A manifold endowed with an affine or metric geometry is said to be geodesically complete 
if all geodesies emanating from any point can be extended to infinite values of the afiine 
parameters in both directions. In the case of a manifold with a positive definite metric it can 
be shown Q that geodesic completeness and metric completeness are equivalent. Focusing 
on a > 2 which is the most likely values of it, the line element Eq. (^) is a vacuum solution 
abstracted away from any source for all values of r. We also realize that 

9 = (1 



r + aM' 

is always positive for all r > 0, so "t" is the invariant world time. The hypersurfaces 
t = const, are spacelike with a positive definite metric 

da^ = (1 '^J^Y^dr^ + (r + aMf{de^ + sin^ ed(P^) (10) 

which is a distance function. Since every Cauchy sequence with respect to this distance 
function converges to a point in the manifold, it is metrically complete. 

Next we show that staticness is not a necessary initial condition if we merely take goe 
independent of time. Let's assume that the line element has the following desired form 

ds^ = B{r,t)dt^ - A{r,t)dr^ - D{r){de'^ + sin^ ^rf^^) (11) 

Taking ds^ = —g^^dx^dx'', the nonvanishing components of the metric connection are 

^, B ^, B' ^, A 

« ~ 25 ' *^ ~ 2B ' ''~2B' 

r-— r-— V"^ - — V - -— r--— sin2/9 

" ~ 2A ' "■ ~ 2A ' 2A ' ^ee - ^A' <t"l' ~ 2A ' 

D' 

r^e = ^ , rj0 = -sin^cos^ , 



2D ' 



rL = cot^. (12) 



where ( ' ) and ( ' ) denote derivatives with respect to r and t respectively. The nonvanishing 
components of the Ricci tensor are 



B" B' A' B' B'D' A A A B. 

" ~ ~2A ^ ilM ^W~ 2AD ^2A~ 4A^A ^ B' 
AD' 



n 



tr 



2AD 



_B" D" D' D' A' B' A' B' A A A B 

''~2B^^~ 2D^^ + A ^ " 4BM ^W~2B^ IbM + W 
D'B' A\ D" 



(13) 
(14) 
(15) 
(16) 



The Einstein field equations with A = for vacuum give 7^^,^ = 0, so from setting 
TZtr = 0, we have 

i = (17) 

and with a reparametrization in the hne element we may also have 

B = (18) 

This means that with the choice Eq. (|TT|) as line element, the solution is necessarily static. 

By considering the bending of light, we try to find an upper bound for a. The deflec- 
tion of light in a gravitational field predicted by the general theory of relativity witnessed 
experimental verification in 1919. Since then, there has been much studies on the gravita- 
tional deflection of light by the Sun and gravitational lensing (GL). Under the great vision 
of Zwicky |]^, observation of a QSO showed the first example of the GL phenomena |^, and 
thereafter it has become the most important tool for probing the universe. It is believed 
that GL can give valuable important information on important questions, such as masses 
of galaxies and clusters of galaxies, the existence of massive exotic objects, determination 
of cosmological parameters and can be also used to test the alternative theories of grav- 
itation 0. The gravitational deflection of light has now been measured more accurately 
at radio wavelenghts with using Very-Long-Baseline Interferometry (VLBI), than at visible 
wavelenghts with available optical techniques. 

Appealing to the spherically symmetric nature of the metric, we consider the geodesies, 
without lose of generality, on the equatorial plane {9 = |). Following Weinberg [^, we get 
the equation for the photon trajectories as 



(r) - 0c 



Air] 



nl/2r 



D r) 



D(r) Birr. 



Dirr,) Br) 



-1 



-1/2 



dr 



The Einstein deflection angle is given by 



A0 



TC 



In Eq. 



we have 
B = (l- 



2M 



r + aM 



) 



A 



2M 



r + aM' 



D 



(r + aMf 



(19) 



(20) 



(21) 



which make the above integral well defined for Tq > (3 — a)M. Since the radial coordinate 
is always non-negative, a merely take positive values. Different values of Tq and a, yield 
different expansions for B and D, so we get different expressions for deflection angle. Our 
investigation is on very small and sufficiently large values of a, since the rest of the interval 
brings no more to us. We give the details of calculating in Appendix A and here use the 
results. For a < 1, the Einstein defiection angle (up to the second order) is as Eq. (|A10|) 



Ax + Ax 



ISvr 



16 



a] 



+ 



(22) 



where 



MG 



X 



and the resteriction imposed by the integral singularity is 

1 
< X < 



3 — a 



thus the closest approach is Tq ^ {3 — a)MG. Switching off a in this solution one recovers the 
well-known Schwarzschild solution, which has extensively examined by many authors (see 
||TT| and refrences therein). However small values is qualitatively similar (but quantitatively 
differnt) to the Schwarzschild case. Further calculation shows that Eq. (p2D is also valid for 
a > 1 (see Eq. ( |A23| )), but it may not contain the closest approach. Therefore, for the weak 
field limit (when r^ is much larger than MG) and all possible values of a, we rewrite the 
equation as 



A0 = A(f)fo 



1 + 



MG /157r 



To 



16 



[l + a] 



+ ■ 



(23) 



where 



Mfo 



AMG 



is the first order defiection angle. The results of VLBl observations of extragalactic radio 
sources, show radio-wave defiection by the Sun ||T2| as 



A0 ^ A0jo(O.9998 ± 0.0008) 



(24) 



Since the order of magnitude of MG/vo for the Sun is 10 ^, Eq. (pHf) and Eq. (|2^ ) make an 
upper bound for a as 



< a < 10' 



(25) 



For sufficiently large values of Q;(actuallya > 3) we may obtain another expression for 
defiection angle (up to the second orther) in the following form (see Eq. ( [A25| )) 



A(/) = By + %^ 



IStt 



-5(a-l) 



+ 



(26) 



where 



and the range of vahdity is 



a{a - 2)MG 



<y < 



a 



thus the closest approach is Tq ^ 0. Eq. (22) and Eq. (p6|) which contain the closest ap- 
proach, can be used for testing the general theory of relativity in a strong gravitational field. 
Although, no test for the theory is known in this region, but there is an open room for such 
investigations. Several possible observational candidates have been proposed to test the 
Einstein theory of relativity in the vicinity of a compact massive object. One of the current 
topics is the study of point source lensing in the strong gravitational field region when the 
deflection angle can be very large ||13|. Our calculation confirms that deflection angle may 
take any small as well as large values depending on a and it would provide a good key for 
the graviational lensing studies. 

Consequently bending of light phenomena, both in the weak field and strong field limits, 
restrict all non- negative values of a as Eq. (|25|) up to 10'^. 



Now we would like to use observational data of the precession of perihelia measurements, 
in order to find a better bound for a. Following Weinberg |T^, for a test particle moves on 
a timelike geodesic in the 6 = 7i/2 plane, the angle swept is given by 



r — r_ 



dr 



^1/2 ( 



D(r) 



D_{B-\r) - Bz^) - D+{B-\r) - B:^^) 
D^D^B^' ~ BZ') 



D(r 



1 1-1/2 

(27) 



where 

D± = D{r±) k 5± = 5(r±) 
The orbit precesses in each revolution by an angle 

A0 = 2|0(r+) -0(r_)| -2n 



(28) 



By using the metric components presented in Eq. (|2l|) , we have gotten the expression for the 
precession per revolution (up to the second order) which the details are given in Appendix 
B. So we may write Eq. (|B14| ) in here as 



A0 = A0y, 






+ ■ 



(29) 



where "L" and "e" are the semilatus rectum and eccentricity respectively, and A(/)/o is the 
well-known first order approximation 



A0/O 



QtiMG 



Fortunately by development of Long-Baseline radio Interferometry and analysis of Radar 



Ranging Data, there are accurate measurements of precession which typically show |15 



A0^A0/o(l.OO3± 0.005) (30) 

Thus matching the theory with observation, with using typical values of ^^ and e, we get 
an upper bound for a as 

< a < 10^ (31) 

We conclude that measurements of these two tests of the general theorey of relativity 
in the weak field limit according to Eq. ( [25| ) and Eq. (pT]), restrict the allowed values of a 
to 10^. Though the observation of the GL phenomena is a difficult task, it would support 
our presented metric components role in a strong gravitational field and would also give an 
accurate bound for a. 

By calculating the Riemann tensor scalar invariant, we receive useful information about 
the existence of singularities. For the line element Eq. (^) in which a 7^ 0, it is equal to 

-oa -o bed _ 48M2 _ 48M2 

^ '^'^- - ^D^ - (r + aMf ^^^^ 

As it is evident, the presence of a makes the scalar finite in the whole range of r, meaning 
that the solutions are free of any intrinsic singularity. Meanwhile there may be a coordinate 
singularity at r = [2 — a)M according to the restriction imposed by Eq. (|25|) . 

We would like to mention two points about this work. One is that, sometimes in liter- 
atures for discussing this problem, coordinate r is defined so that the area of the surfaces 
r = const, to be Airr'^ [|16|. This generally is not the case, since before fixing the metric 
there is no possibility of speaking the distance, then in the same way, there is no possibility 
of speaking the area. We take r = (x^ + y^ + z'^Y^'^ where (x, y, z) are usual cartesian space 
coordinates. 

The other point which should be taken with caution is, while at first r was taken as a 
space coordinate with the range (0, 00) and the particle was located at r = 0, at end we 
come to the conclusion that r is merely a space coordinate in the interval ((2 — a)M., 00). 
For the rest of the interval (0, (2 — a)M), it is standing as a time coordinate. This contra- 
diction or at least ambiguity raises the question that While the location of the particle is not 
well-defined, how may we speak of the value of D at this point? This ambiguity should be 
found a satisfactory explanation. 



III. CASE A / , a = 

Recently for vaccum sphericaly symmetric space, a non-static solution of Einstein field 
equations with cosmological constant in a proper FRW type coordinate system has been 
cosidered [^]. The result shows a singularity at the origin which the intrinsic nature of 

8 



it may be checked by calculating the invariant Gaussian curvature or the Riemann tensor 
scalar invariant. In this section we choose tensor analysis to find the answer. So we start 
with the following form of the line element 

ds'^ = B{r, t)dt^ - R\t) A{r, t)dr'^ + r'^i.dO'^ + sin^ edct)"^) (33) 

where R{t) is the scale factor and satisfies the Friedmann equation 

R^ R ^ 

and A is the cosmo logical constant. The cosmological time t, which determines the evolution 
of the universe is measured by the proper time of clocks fixed on the geodesically moving 
galaxies comprising the universe, synchronized such that the time t = corresponds to the 
beginning of the universe at the big bang. 

According to Ref. |]3] the nonvanishing components of the metric connection are 



r* - 
^ tt — 


B 
2B ' 


r* - 

^ tr — 


B' 
2B ' 


r* = 

rr 


RRA R^A ^, 


RRr^ 
B 


F* - 


B ^^"^ 


v — 
^ tt — 


B' 


-rr _ 
^ tr ~ 


R 

R + 


A 
2A' 


A' 


r 




-^sin^^, 


2AR^ 


J- te — 


R 
R' 


p9 - 


1 
r 




rJ^ = -sin^cos^, 








^ t<t> — 


R 
R' 


'p4> _ 

'- rcf> - 


1 
r 




rt = cot^. 









(35) 



The nonvanishing components of the curvature tensor then become 

^ _ AR^ AR^ B" 3ARR _ B" A 3AR {Rf 



B 2B 2B 2B ' 2AR^ 2A 2AR R^ 

, A'r r^RR A R ^ 

/<■ ete - -^ H '^^2A "^ ^ ^''^ ~ ^*^ ' 

T^ M(h = T^ ete sin 6 , TZ^ ^rA = TZ d>t<h ■, 



B R B R^ 

2rAR^ '^R^W^R' ' 
A' RARR R^A^ 



^ tot - o^/.D2 + -b(w^ ^ ^) ' ^"^ t<pt - 7^^ tet 



/C rSr — „ . + p (, „ + j , /v r<j>r — 'V. rOr , 

1 r^/?2 

n' ^e^ = TZ't' e^esin^^ , 7^'/' ^^^ = 1 _ + . (35) 

A B 

Since the metric components have no cross terms and furthermore the curvature tensors 
have some algebraic properties, the Riemann tensor scalar invariant takes a simpler form as 

7^'^ bcdTla ^'^'^ = 2^2 ig^'^T^" m-m)' (37) 



Using Eq. (^6]) gives 

^j~) a ij-) bed 
K- bcdK-a 



B" A _ 3AR i?2 Y 

2ARR2 ~ ^Xb ^ 2ABR ~ 'BR^J 



+ 16 



A' 



+ 



AR 



+ 



R^ \ 



2rA^R^ ' 2ABR ' BR^J 
/ 1 1 R^ Y 

The rr component of the field equation gives 

TZrr + ^9rr = 

and makes the equall terms for the first parentheses as 



B" 



A 



3AR 



i?2 



2ABR^ 2AB 2ABR BR^ 
The 66 component gives 

TZeg + Aggg = 

which this leads to 



1 1 R^ \ 



A' AR i?2 



1 



1 i?2 3/^2 

+ 



2rA^R^ 2ABR BR^ 2\r'^R'^ r'^AR^ BR^ R^ 

By solving field equations we also have [§] 

2M 



l-\ + ^r^R\A-l) 



rR 



or 



B = A-^ = - 
2 



ri? 3 ^ 3 ^ 



2M 

7r 



3 ^ 



Inserting 



(13) and (|T]) in 



"^^ bcd'^a 



bed 



gives 

48M2 



24,2 
\ A 

r^R^it) 9 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



This evidently exhibits the existence of an intrinsic singularity at the origin. Removing this 
deficiency, leads us to the most general form of the solutions which comes next. As it is 



expected, Eq. ( ^31) with R{t) = 1 gives the result of the static case |T7 



We end this section with emphasizing on some aspects of Eq. (42) and Eq. (|43|) . Firstly, 
the existence of a nonzero cosmological constant regardless of its actual value, is sufficient 
to prevent from Schwarzschild singularity ai r k. 2M. Recent observations of type la 
supernovae indicating a universal expansion, put forward the possible existence of a small 



10 



positive cosmological constant [Q. These evidences persuade us that in a A-dominated 
universe, we would have no troubles in describing the whole space. 

Secondly, since there is no singularity for r > 0, then there is no ambiguity in defining 
coordinate r, which mentioned at the end of Sec. 11. "r" is a space coordinate in the whole 
interval (0, cxd) and we may speak of the value of D at r = without any problem. 

Thirdly, a coordinate transformation |^, transforms the above non-static metric to the 
well-known Schwarzschild-de Sitter static metric Eq. (j^). This metric has some deficiencies 
which we would like to discuss briefiy. Despite of an intrinsic singularity at the origin which is 



similar to our case, it has a coordinate type singularity at r ^ Jj^ in a A-dominated universe. 
Though the presence of cosmological constsnt removes the Schwarzschild singularity from 
our metric, but the transformed form of it, i.e. Schwarzschild-de Sitter metric, has again the 
problem of exchanging the meaning of space and time. On the other hand when M = 0, the 
assumed FRW background could not be revisited but the attached homogenity and isotropy 
would resist. And more importantly, this metric shows a redshift-magnitude relation which 
contradicts the observational data [^]. 

Therefore it is adequate to discard the Schwarzschild-de Sitter metric, in favor of our 
presented metric, as a proper frame of reference in the presence of A, since all of the men- 
tioned deficiencies do not exist in it. The metric asymptotically approaches to the non-static 
de Sitter metric which is appropriate for a A-dominated universe. Furthermore as we show 
next, the presented metric has the suitability even to remove the intrinsic singularity at the 
origin. 



IV. CASE A / , a / 

Since it turns out that there is an intrinsic singularity with the choice a = 0, we would like 
to solve the problem by using the most general form of the line element. In this section we 
work out expressions for metric coeficients in the FRW universe and obtain an analytic metric 
everywhere. Therefore we choose the metric for the universe in terms of the coordinates (r, t) 
to be 



ds^ = B{r,t)de - R^{t) 



(44) 



A(r, t)dr'^ + D{r){de^ + sin^ ' 

The nonvanishing components of the metric and metric connection respectively are 

gu = -B, grr = AR\ gee = DR^, g^^ = DR^ sin^ 6 (45) 

B „, B' „, RRA R^A , _ RRD , _ RRD . ^ 

r^<^ = -sin^cos^ , 

rt = cot e. (46) 

11 



r* - 

'- tt — 


2B ' 


r* - 


F* 

rr 


-rr _ 

^ tt — 


B' 


^'^~ R^ 


A 


2AR^ 


2A 


J- te — 


R 
R' 


^^' 2D 




-p't' — 
^ t<i> — 


R 
R' 


r4> 2D 


•) 



where ( ' ) and ( ' ) denote derivatives with respect to r and t respectively. The nonvanishing 
components of the Ricci tensor then become 



ntt 



B 

AR2 



B" B'A' B'. B'D' 



7^ 



2B AB^ A B' 2BD 
3R R 3B A 
'B^~R^2B~A^ 
AD' B'R 



A A A B, 

2A~4A^A^B' 



tr 



2AD BR 

B" D" D' ,D' Ar. B' ,A' B\ 



/\yf"P — .^ I .^ 



2B D 2D'D A' 45 M B' 
R^A\R R 2R 2A B A A A B, 
.R ^ R^li ^ ^ ' 2B' ^ 2A~ 4AM ^ 'B' 



B 



TZrhrh = sin^ 971, 



D'B' A\ D" R^D 



^'^AA^B A'^2A B 



R R 2R A B^ 
'R^R^^^2A~2B' 



(47) 
(48) 

(49) 
(50) 



To solve the vacuum field equations TZ^i, + Ag^i, = 0, we first begin with TZtr = and 
introduce a new variable p to be 



p = R{t)D^I\r 



(51) 



From this equation we obtain 



1 4* R* 

2 ^ A B' 



(52) 



which (*) means differentiation with respect to p. Since R and D' ^ 0, we should have 

A* B* 



(53) 



Integration with respect to p and imposing fiat boundary condition at large distances yields 

AB = 1 (54) 

in agreement with the FRW background. In the next step, let's consider 



B AR^ 







(55) 



which we get 



-D 



/2 



D" 



or we may write 



2AD2i?2 ADR^ 



{D'D-^/^y = 



(56) 



12 



Integration with respect to r gives 

D^/'^ = r + aM or D = {r + aMf 

As before a is a positive constant in the range < a < 10^. 

Finally from the 99 component of the field equation, TZqq + Agfg 
functional form of A{r, t) as follows 

-1 + 1 - p^ - kp\A - 1) - ^p^A* = 



(57) 



0, we obtain the 



or 



d 


- 


1 1 


A d 




p('\ - 




H 


dp 




A'\ 


3 dp 



Integration then yields 



pHa-1) 



p{l'j) + ^p'iA-l) = c 



(58) 



(59) 



where c is constant of integration. Comparing our result with post-Newtonian limit gives 
c = 2M. Finally the result is 



A = B-^ 



,2A 



,2\-l 



2M A 4A 2M A ,. 
.1 tP r + -irP - 1 tP) 

p 6 6 p 6 



or 



B = A-^ 



2M A 



P 



-P 



2^2 



4A 



p2 + (1 _ 



2M A 



P 



-P 



(60) 



As it is evident from the functional form of A, B and D, this metric has no apparent 
singularity and a straightforward calculation gives the Riemann tensor scalar invariant in 
the form 



Tl"^ bcdT^a 



bed 



48M2 24 , 2 
H A^ 



Dm^{t) 9 



(61) 



If we set A = 0, R{t) = 1 in (|6lD, Eq. ( 13^) will be obtained, and furthermore with D = r^, 
Eq. (H) would regain. 

It is remarkable that having a 7^ 0, we would have a well-defined metric in the whole 
space which asymptotically approaches to the non-static de Sitter metric; i. e. the appropri- 
ate metric for a A-dominated universe. 
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V. GEODESIC EQUATIONS 

Our next task is to obtain and solve the geodesic equations of a freely falling material 
particle in a proper FRW type coordinate system. In the static case, the basic equations 
which determine the geodesic structure of the manifold have been fully developed |[r?| , p!8 . 
We solve the problem in the non-static case by using the equations of free fall 

(fx'^ ^,, dx" dx" ^ ,„„, 

Using the nonvanishing components of afiine connection, given by Eq. (^H]) we have 

dH B ,dt,y B'dtdr RRA R^A.Ar,^ RRD AO ,^ 
ds"^ 2B^ds' Bdsds ^ B 2B'^ds' B ^ds' 

d'^r B' dt 2 r,fR A dtdr A' dr 2 D' dO 2 

d?^ 2AR'^^~ds ^ ^R^2A}Ts~ds^2A^Ts' ~ 2A^~ds 

_:^sin2^(^)2 = (64) 

2A ^ds' ^ ' 

d'^e Rdtde D'drde _ . ^/^0n2_q /(..n 

ds^ Rds ds D ds ds ds 

^ + 2- — — + — — — + 2cot^— — = (66) 

ds^ Rds ds D ds ds ds ds 

Since the field is isotropic, we may consider the orbit of our particle to be confined to 
the equatorial plane, that is 6' = |. Then Eq. ( |65D immediately is satisfied and we may 
forget about ^ as a dynamical variable. Introducing the previous variable "p" and using the 
Friedmann equation, make the above equations as 

d^p [a dH A A* A 2A* dt2^A*dp2 pM.2 n (a^\ 

d'^cj) 2dpd(f) , . 

ds'^ p ds ds 

The last equation gives 

Integrating ( [70D with respect to s we get 
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where J is constant of integration. 

In order to derive the exact solution, first rewrite Eq. ( |67D as 

Ik dH , A 



d_ 
ds 



, 1 A ^ ^.dt' 



+ \l^P^ 



3^rfs2+4A3 



A A ^A\,dt^^ J2 
^P + TP-r)(-r) + 



6 A'^ds' 



Ap3 



Wt(^-f-)(S'^ 







Substituting Eq. (H) into Eq. (^) yields 



ds 



1 A ^ ,,dt 



+ 



rfV 



^*^,dp, 



pA^^UA + pA-)(fJ'^ 



or 






A A o ,, (it /A ,(ip' 



which can be integrated to get 

,1 A ^ ^.dt [a dp 

where Ci is constant of integration. Derivative of p with respect to s gives 

dp A dt dr 

ds V 3 (is ds 

if r — >• oo or p — >• oo, then (|75[) reduces to 



dp 

ds 



-P 



1 



(72) 



(73) 



(74) 



(75) 



because for free fall at infinity we have 

ds ds 

Therefore Ci would be fixed by asymptotic behavior as Ci = 1. Eq. ([7^) then becomes 



dt _ 1 
ds 



fp^S 



1 _ An2 4 



(76) 



We have yet another equation, the line element equation, which in terms of p variable 
has the form 



_ :^p'A){^f + 2J'-;pA^^ - aA' - p2(^)2 = 1 



A 3 



.dt . 
'Ts' 



I A dt dp 
3 ds ds 



ds' 



ds 



(77) 
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Rewriting Eq. 
(ft 



as 



A A* dt 2 A* (it dp 



ds^^^ 3^2A'ds 



.JIa 



A 
3"' 



^\r-- 



A ds ds 

dt 2 r, l^ .dt dp 
'ds V 3 ds ds ' '' ds' 

I A dt dp 



ip'M^y-^\rip^T-i^+M'^r+'4 



p- 



^ 2A(d't.2 9 /^ dt dp A(dP\2 
3 ds \ 3 ds ds ds . 



and using Eq. ([77|) , we come to 
dH lk,_ A\,dt 



2 A* dt dp 



ds"^ V 3 A ds A ds ds 



'A(a+|a* + ^A* 



multiplying the above equation by A/f p and inserting the result in Eq. ( |68D we find 



(78) 



(79) 



d^p A* 



ds^ 2^2 



, 1 ^ 2 A\/dt.2 r^ /^ .dtdp .,dp,o 

^A-r^^^cLs^^ir^d-sTs-^^Ts^. 



p^A 



-^p(A + ^A* + —A*) 
3'^^ 2 2p ' 







(80) 



From Eq. (|77D , we notice that the expression in the bracket is just equal to 1 + ^, so we 
obtain a second order differential equation for p as 



'1 + -T 



P 



ds^ 2^2 
which can be written 

d'^p 1 d 



or 



ds'' 

d 
ds 



+ 7:^ 



4L_Ap(A+^A* + — A*) 
P'^A 3^^ 2 2p ' 



i\-TP'm + ^) 



2dpVA 3 



P^ 



:|y^+(^-?r^)(i-H^ 



1 A 
I 3' 



J_2 



Integration then yields 



:'£r + (i - f .M)(i . i!; 



C2 



Asymptotic behavior shows that C2 = 1, and thereby 



,dp 



1 A 



J^ 



Eq. ( ^91) together with J = makes Eq. (^21) as 



(81) 



(82) 



(83) 



(84) 



(85) 
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<Pp d 
ds"^ dp 



-M A 



P 



6 



P 



(86) 



Since the potential field $ is defined as 



then 



$ 



ds"^ dp 



-M A 



P 



r 



(87) 



which is equivalent to the Schwarzschild-de Sitter potential [|T9|. Although the potential at 
r = i.e. p = aMR{t) is very large, but it is finite. It seems likely that the potential field 
of massive stars show this behavior and help us to find a physical mechainsm for fixing a. 
It will be a great success if observing extra high energy phenomena in AGN's and cosmic 
rays, provide a lower limit for a. 



REMARKS 

A space-time is said to be spherically symmetric, if it admits the group SO (3) as a group 
of isometrics, with the group orbits spacelike two surfaces. A coordinate transformation like 
r —>■ r' = r + aM which translates the center of symmetry and thereby breakdowns spher- 
ical symmetry does not belong to SO (3). We have shown that the solutions of spherically 
symmetric vacuum Einstein field equations are necessarily neither static nor incomplete. 



APPENDIX A: DEFLECTION OF LIGHT 



We would like to derive the deflection of light formula which introduce in Sec. II. The 
orbit is described by Eq. (p!9|), that is 



A{r) 



D(r) 



1/2 



D(r) B{r„ 



D(r^ B(r] 



-1/2 



dr 



(Al) 



Using for A{r), B{r) and D{r) the Eq. (|2TD, we notice that the integral is defined for 
To > (3 — a)M. The first term of the second square root is written as 



D{r) Bjro) 
Dir^ B(r) 



aM/r 



1 + aM/vo 



- (2 - a)Mlr 
(2 - a)Mlr, 



(A2) 



We want to derive deflection angle up to the second order but the expansion explicity de- 
pends on a and M/to. Since there is no common valid range of expansion for the involved 
parameters, we must treat the necessary and sufficient cases seperately; i.e. a < 1 and 
a > 1. 
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We consider first a < 1. Because of tlie integral singularity condition, we always have 
To > (2 — a)M and therefore with a simple expansion we get 



D{r) B{ro 
D{r,) B{r) 



l + (2(l + a)M 



+ 



(2-aYM^ n 1 



r rr 



6aM^ 



r \r ro 

The argument of the second square root in Eq. (|A1|) is then 



1 1^2 



r To 



D{r) B{ro, 
D(r^ B(r) 



-1 



-) -1 

r, 



{2{l + a)M-3a'^Myro)r 

ro{r + To) 
{2-afM^ 6aM'^{r-ro) 



Toir + ro) rl{r + ro) 



so it gives 



Dir^ Bir) 



-1/2 



-) -1 

r, 



-1/2 



2/\/r2 



(2(l + a)M-3a2MVro)r (2 - a)2M' 



2ro(r + To) ' 2ro(r + Tq 

3aM2(r-ro) 3(1 + a)2MV 



'^ol^ + '^o) 2r2(r + ro)^ 



+ ■■ 



The argument of the first square root in Eq. ( |A1| ) and its expansion is 

1-1/2 

(1 + aM/r){l - (2 - a)M/r) 



M(r)l 


1/2 _ 1 


[D{r)\ 


r 




1 




r 



(l-a)M (2a2_4a + 3)M2 

Inserting these two equations into Eq. (|A1|) gives 

dr 



{ro) - 0c 



. r[(f-)2 - l]i/2 



(l-a)M (2a2-4a + 3)M2 
r 2r'^ 



+ 



(2(1 + a) M-3a'^M^/ro)r {2-afM'^ 
2ro{r + To) 2ro{r + r^) 

3aM'^(r-ro) 3(1 + a)2MV 



'^o('^ + '^o) 2r2(r + ro)^ 



+ 



The integral is elementary, and gives 



TT 2M 2M^\lhTx , 

77 + + —IT- ^TT - 1 + a 



°° 2 ' To ' r2 
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+ ■■ 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



The deflection of the orbit from a straight line is 



TT 



A0 = 2|0(ro)-0o,, 
Hence to second order in M/to-, the deflection angle is as follows 



A0 



AM AM^ 



IStt 



-(l + a) 



+ 



(A9) 



(AlO) 



Next we consider a > 1. Since in this case we may also have Tq < |2 — q;|M (e.g. for 
a > 3 the integral singularity condition reduces to Tq > 0), so in order to have a valid 
expansion we define the following parameters 



aM 



& 



(2-a)M 



1 - (2 - a)M/ro 



1 + aM/vo 
and with a > 1, we always have 

a/ro < 1 & b/ro < 1 
So we rewrite Eq. (|A2|) and expand it as 



(All) 



D{r) B{r,) 
Dirr) Bir) 



■ ^ ■ 


2r 


To. 




- J. ■ 


2 


To- 





, 1 1 

l + a( 

r To 



l3 



1-6 



r To 



T-l 



1 1 



1 + (3a + 6) ( ) + (3a2 + 3a6 + 6^) ( - - — ) + ■ ■ 



The argument of the second square root in Eq. ( |A1| ) is then 



Dir) Biro) 



1 



D{ro) B{r] 
so it gives 

Dir) B{ro) 



(3a + b)r {3a'^ + 3ab + b^){r- r^) 



ro{r + ro) 



rl{r + ro) 



+ 



(A12) 



(A13) 



Dir^ Bir) 



-1 



-1/2 



-1/2 



(3a + 6)r 3(3a + fe)^^ 
2ro(r + ro) 8r2(r + ro)^ 
(3a2 + 3a6 + 62)(r-ro) 
2r2(r + To) 



The argument of the first square root in Eq. (|A1|) is 
■A(r)lV2 1 



Dir) 



[l + aM/r){l-{2-a)M/r) 



-1/2 



In terms of these new parameters, we have the following equalities 

1 + aM/r = il-a/ro)'^ 
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l + a( 

r rn 



(A14) 



(A15) 



(A16) 



and 



1 - (2 - a)Mlr = {l + hjroY^ 






(A17) 



Inserting these results into Eq. ( [A15| ) and expanding it up to second order yields 



M(r)l 
\d[t)\ 


1/2 


1 ■ 

r . 

X 



(l_a/r„)(l+Vro) 

1 _ (^-^) (\ _ M 
2 Vr To) 



1/2 



(3a2 - 2a6 + 362) .^ ^.2 
8 Vr r^, 



(A18) 



Substituting Eq. (|A14| ) and Eq. ( |A18| ) into Eq. (|A1| ) gives an elementary integral as 



[To) - 0c 



dr 



1 



nl/2 



X 



X 



{l-a/ro){l + h/ro) 



a-h) (I 1\ (3a2-2a6 + 362) /I 1 



r r^ 



r r^ 



+ 



(3a + 6)r 3(3a + 6) V {"ia^ + ^ah + }p){r - To) 



2ro{r + ro) 8r'^{r + VoY 



2rl{r + ro) 



which can be easily integrated 



_ vr {a + h) (a + 6)(a + 36) 157r(a + 6)^ 
0(r„) - 0^ _ - + — - ^-^ + ^^-^ + 

This gives the deflection angle in terms of the defined parameters as 



(A19) 



(A20) 



A0 



2{a + h) (a + 6)(a + 36) 157r(a + 6) 



+ 



16r^ 



+ 



(A21) 



For expanding the parameters, we must notice the region in which the expansion is valid. 
For To > aM (and consequently r^ > |2 — q;|M) the expansions to second order in M/to are 



a/ro 



aM a^M 



2/\#2 



+ ■ 



Tn 






2 - a)M 2 - afM^ 
h To = '— + ^ { + ■ • • 

Hence the Einstein deflection angle up to the second order in this region becomes 



(A22) 



A0 



AM AM 

+ ■ 



2 r 



To 



ISvr 



(1 + a) 



+ 



(A23) 



The other region we are interested in, is Tq < |2 — a|M (Actually this region exists if 
a > 2.5) which gives the second order expansions as 
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h/ro 



1- 



r„ 



+ 



TZ, 



-1 + 



To 



+ 



+ ■ 



(a - 2)M (a - 2)2M2 
Therefore the Einstein deflection angle up to the second order in this region is 

8ro 4r2 flSvr 



(A24) 



^'^ a(a-2)M^a2(^_2)2M2 
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-5(a-l) 



+ 



(A25) 



APPENDIX B: PRECESSION OF PERIHELIA 

In order to derive the Precession of Perihelia expression which is used in Sec. II, we 
consider a test particle bound in an orbit around a massive object. The angle swept is given 
by Eq. (Ig^), that is 



r_ 



D{r_) {B-\r) - B-\r_)) - D{r+) {B-\r) - B-\r+)) 1 



X 



Ml/2(, 



D(r] 



D{r+)D{r_) iB-^{r+) - fi-i(r_)) 



dr 



D(r) 



-1/2 



and the orbit precesses in each revolution by an angle 

A0 = 2|0(r+) -0(r_)| - 27r 



(Bl) 



(B2) 



Following Weinberg, we make the argument of the first square root in Eq. (pT| ) a quadratic 
function of 1/D{r) which vanishes at D{r±), so 

D{r_) (5-i(r) - B-^{r_)) - D{r+) {B-\r) - B-^{r+)) 
D{r+)D{r^) iB-^{r+) - 5-i(r_)) 

1 " ( '^ 1 \ / 1 1 A 



C 



The constant C can be determined by letting r ^ oo 

D{r+) (1 - B-\r+)) - Djr,) (1 - B-\r.)) 
D^/^{r+)Dy^{r_) {B-\r+) - B-\r_)) 

and its expansion up to second order becomes 

C = 1 - 2M ( — + — 1 + 2aM2 (\ + \\+... 



(B3) 



(B4) 



(B5) 
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Using Eq. (0) in Eq. (^ gives then 



r — r_ 



^-1/2 r 

Jr- 



dr 



AV2( 



D(r] 



)(- 



^Dl/2(^_) £)l/2^^-)^V^l/2(^^ JJl/ 



Intoducing a new variable ip 
1 1/1 



+ 



makes the integral as 



sin-?/; 



ir) - 0(r_) = C"^/^ / (i'?/' 

-7r/2 



1-M 



-M 



1 



1 1 



sin-?/' 



-1/2 



Expanding the square root to second order in M/D^/^(r), and integrating yields 



(r+) - </)(r_) = C-^/2 



7r+ 



ttM 



1 1 

+ 



2 V^^/^K) D^l^{r^)) 

1-2/-. ^ \ 2 



+ 



37rM^ 



8 \D^I\r 
SttM^ / 1 



1 1 

r + 



Dl/2(^_] 



1 



16 \D^I\r^) D^/^ 
Using Eq. ( p5|) in the above equation and expanding it gives 



+ ■■ 



V r_i_ r_ 



- aM^ — 



X 



ttM / 1 1 ' 

7r + ^- — + — 
2 V r 1 r_ , 



2 \^r_|_ r^J 

naA'P /I 1 ' 

2~ld!'^^, 



rl r: 



+ 



SttM^ /I 1 \ ^ SttM^ /I 1 ' ^ 



^r+ r^J 16 \r+ r_y 



So with Eq. ([B^, we get the precession per revolution as 



, / 1 1 \ 197rM2 /I 1 ' 
SttM — + — + ^ — + — 

-SvraM^ [-^ + -^] + • • • 



SttM^ / 1 



r I r_ 



The elements of planetary orbits are a and e which defined by 



-1/2 



(B6) 



(B7) 



fB8) 



(B9) 



(BIO) 



(Bll) 
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r±= (1 ± e)a 
L=(l-e2)a (B12) 



where 

1 1/1 1 



These result the desired form of precession per revolution up to the second order as 



A0 



GttMG 



2 



MG / 19 e , ., 

l + ^(- + --2a(l + e- 



+ ■■• (B14) 
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